Proof that the Hydrogen-antihydrogen Molecule is Unstable. 
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In the framework of nonrelativistic quantum mechanics we derive a necessary condition 
for four Coulomb charges (rri±, , ni^, ), where all masses are assumed finite, to form 
the stable system. The obtained stability condition is physical and is expressed through the 
required minimal ratio of Jacobi masses. In particular this provides the rigorous proof that 
the hydrogen-antihydrogen molecule is unstable. This is the first result of this sort for four 



particles. 
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Recent success in the production of trapped antihydrogen atoms |l| has renewed interest in the 
g ; interaction of matter with antimatter and especially in the hydrogen-antihydrogen system (H-H). It 

has long been conjectured that with pure Coulomb forces no bound state of hydrogen-antihydrogen 
exists. The numerical evidence supports this conjecture yet there is a lack of rigorous proof 



as remarked by some authors 
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J. Our aim in this Letter is: (i) to supply such proof under 
assumption that only Coulomb forces act between the constituents, (ii) to provide insight into the 
screening effect within the system of four charged particles. 

The instability of H-H is explained by the screening effect just like the instability of the muonic 



hydrogen ion (H/x ). In the system (H/i ) the heavy muon gets so tight to the nucleus that screens 



the positive charge and the electron "sees" a tightly bound neutral combination pn~ and departs 
from it making the whole system unbound. In [5( we have proved that the screening effect in the 
system of three charged particles is expressed through some critical ratio of Jacobi masses. From 
the physical point of view there are two orbits in this system to consider, namely one orbit within 
the pair of particles (the pair that sets up the dissociation threshold) and the orbit of the third 
particle in the field of this pair with respect to the pair's center of mass. Inverse Jacobi masses 
are proportional to the Bohr radii of these orbits. If, say, the orbit of one negative particle is 
outdistanced then the attraction from the positive charge is screened by the other negative particle 
and the system becomes unbound. 
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The system of four unit charges (mt , , rri^ , m J ) can be unstable only against dissociation 
into two neutral pairs. Indeed, if the lowest dissociation threshold would be dissociation into 
one particle and the bound cluster of three particles, then these two objects would have opposite 
charges and the long-tailed Coulomb attraction between them would guarantee the existence of a 
bound state below the dissociation threshold. (Just the same argument explains why atoms are 
stable). This suggests that we have to consider three orbits, two inner orbits of the neutral pairs 
and the third orbit of the relative motion of these pairs. The Jacobi masses for the neutral pairs 
are \i x = mfm^/imf + m^~) and \i y = /(m^ +mj). The Jacobi mass corresponding to the 

relative motion of these two pairs is [jlr = (mf + m^~)(m^~ + m 4 )/(m| + + + mj). Pay 
attention that we order the particles so that among two possible rearrangements into neutral pairs 
the lowest energy threshold corresponds to the dissociation into (m^m^) + (m^~mj) and the pairs 
are ordered so that \i x > fj, y . 

The Jacobi masses are in fact not independent, it is easy to check that (Ir > 4^ y always holds 
if the particles are ordered as above. Let us consider the screening effect within the system of four 
particles keeping in mind that the Bohr radii of the orbits are inverse proportional to the Jacobi 
masses. The first possibility is fi y <C ^r and fiR ~ \x x . This would mean that three particles form 
a heavy cluster and one particle is outdistanced from it. Here everything depends on whether this 
heavy three-body cluster has a bound state. If it does then the whole system is stable because the 
cluster and the particle have opposite charges. Hence the whole situation reduces to the question 
whether there is any screening in the three-body system. For ur ~ \x x there is no apparent 
screening as follows from 0,0. For example, Bressanini et. al. P] have collected the convincing 
evidence that the system (mf, 1~ ,mf, 1~) is stable for any and m^. The three-body system 
(m^, 1~ , 1~) is always stable and if <C 1 then we run into the situation where n y <C ^r and 
fJ-R ~ and still the whole system is stable. This means that ^ y fiR is not sufficient for the 
screening effect to take over. 

Let us consider the possibility when ur <C which would mean that the pair (mf, m^) has a 
very short inner orbit and other particles are outdistanced from it. In this case it is right to expect 
screening because other charged particles would "see" the tightly bound pair (m^ , ) as neutral 
and the system would fall apart. Our aim in this Letter is to present the rigorous and analytic 
proof of this screening effect, namely 




Instability, 



(1) 



where under Instability we mean the absence of a bound state below the dissociation thresholds. 
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Eq. (J2) manifests the screening effect for four particles. From Eq. (^Q) it easily follows that the 
hydrogen-antihydrogen molecule has no bound state and must decay into protonium and positro- 
nium. Muonic molecules pn~e + e~ and fi + pe + e~ are unstable as well. 

The proof of Eq. (^Q) is along the same line as in 5] (the basic idea of the proof is due to Thirring 
3]). Before we proceed with the proof let us introduce the notations. Let qt, r-j 6 M 3 denote charges 
and position vectors of the particles i = 1,2, 3, 4. We shall work in the system of units where H = 1. 
We put qi t z = +1, and #2,4 = — 1, and the interactions between the particles are = q%qk/\^i — *k\ 
(remember how the particles are ordered). The stability problem with Coulomb interactions is 
invariant with respect to scaling all masses f|, so we can put [i x = 2. By the end we shall rescale 
the masses back. 

Consider the system of four charged particles which is stable for fiR < 3/8 (this is weaker than 
in Eq. l[T])). We separate the center of mass motion in the Jacobi frame ,§t] putting x = T2 — ri, 
y = r-4 — r-3, R = —ax + r-3 — ri + by, where a = 7712/(7711 + 771,2) and b = 777,4/(777-3 + 777,4) are 
the mass parameters invariant with respect to mass scaling. With Jacobi momenta defined as 
Px,y,R = —i^x,y,R the Hamiltonian of the system takes the form 

H = hi 2 + h u + ^- + W, (2) 
2/4R 

where 

W = Via + V U + V 23 + V24, (3) 
and hi2 = P^/4 — 1/x, and /134 = Py/(2fi y ) — 1/y are the Hamiltonians of the pairs (1,2) and (3,4) 



(notation x is used instead of |x|). The ground state wave function of /t, 12 is 4>q = Y / 8/7rexp(— 2x) 
so that hi2<t>o = —(fro- By the particle ordering the energy threshold corresponding to dissociation 
into two neutral pairs is E t h = —1-/^/2, which is the sum of the binding energies of the pairs (1,2) 
and (3,4). Following 5] we shall cut off the positive part of W by introducing W- = (\W\ — W)/2 
and W + = (\W\ + W)/2 which results in the decomposition W = W + — W-, where W± > 0. 
Instead of H we shall consider the Hamiltonian 

H = h 12 + h 34 + ^-W„. (4) 

n 

(The operator H is self-adjoint on the same domain as H, see ,5]). We shall assume that H is 
stable, i.e. H has a bound state $ with the energy E < Eth- Because H < H we conclude that 
(^1^1^) < £ , i / l ||^I , || 2 . Before we use this inequality let us introduce a projection operator P , which 
acts on any /(x, y,R) as 

P f = M%) J dxVo(xO/(x',y,R), (5) 
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and put 7] = Pq^ and £ = (1 — Pq)^, where obviously 77_L£ and = rj + £. We shall assume that 
1 1 £|| ^ (later we shall get rid of this assumption), then we are free to choose such normalization 
of \F that ||£|| = 1. Now let us rewrite the inequality (\I'|1I|\E') < E^H^H 2 decomposing ^ into 
^ = r/ + £. 

(n\hu\ri) + Hp- - w-\n) - (r/lWLio - &\w-\n) 

2^r 

+<*IM0 + <f M*) - (t\&. - W .\0 < -1 - (Hy/2)(\\ V f + 1), (6) 

where we have used that the terms like (jy|Py|£) cancel because Pq commutes with the operators 
Py,p|j and 1/y. Indeed, in this case for example (r/|Py|£) = {v\PopI\0 = (v\PyPo\Q = 0- 

We are going to rewrite Eq. © using lower bounds for some of its terms. From the Hydrogen 
ground state and by the variational principle for the terms in Eq. © the following inequalities 
hold (r]\h34\rj) > — (/J-y/2)\\r}\\ 2 and (£|/i34|£) > —fJ, y /2. Introducing two non-negative constants 



a = \J (rj\W-\rj) and (3 = y (£|W_|£) we get by virtue of the Schwarz inequality |(^|W_|r/)| < a/3. 
It remains to figure out the bound for the term {t,\h\ 2 \£). From the bound spectrum of the Hydrogen 
atom we have we have [?J h\ 2 > — Pq — 1/4(1 — Po). (Indeed, Pq projects on the ground state of 
hi2 which has the energy —1 and the energy of all other excited states is greater or equal to —1/4 
which is the energy of the second excited state). Hence for the first term in Eq. © we get the 
bound (£|/ii2|£) > —1/4. Substituting this into Eq. (jSJ) leaves us with the main inequality 

{rjl ^ ~ WAV) ~ 2a/? + (Cl 2^ ~ W ~ 10 < "I (7) 

We shall focus on the third term on the lhs of Eq. Ij7]). 
First, let us prove that the inequality 

+ AV 14 + AV 23 > -2A 2 m (8) 

2^R 

holds in the operator sense for any A > 0. Indeed, the interactions have the form V14 = — 1/|R — z±\ 
and V23 = — 1/|R — z 2 \, where the vectors z± = —ax — (1 — b)y and z 2 = (1 — a)x + by play the 
role of parameters in Eq. (JSJ). According to the result of Lieb and Simon [l| the energy of the 
Hamiltonian on the lhs of Eq. (JHJ) is monotonically increasing with |zi — z 2 \. Hence the minimum 
energy is attained when z\ = z 2 = 0, which gives us Eq. (JSJ). From Eq. (JHJ) using the obvious 
inequality — W- > V14 + V23 we find that for any A > and x( x > y> R-) 

( X \J?L - AW.\x) > ^A^nWxf (9) 

2flR 
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holds. With the help of Eq. @ we obtain the following chain of inequalities 



<^-^_|e> = max 
Z[l R A>-1 



' > (io) 



(£1^-^+1)^-10 + A/3 2 



max [-2(A + I) 2 m + A/3 2 ] = A- - /? 2 , (11) 

A>-1 

where in Eq. (jl(J|) we have added and subtracted the term A/3 2 = A(£|W_|£). Substituting Eq. I|1U|)- 
((TT|) into Eq. (JJJ) leaves us with the inequality 

(r? |^i _ p^fo) + A. _ ^ _ 2a/3 < _3 (12) 
The following inequality always holds 



To see that Eq. (|13|> is true it suffices to take all terms to the lhs and minimize over a. Substituting 
Eq. (|13|) into Eq. (|12j) and using a 2 = (r/|W_|7/) allows us to formulate the stability condition 

(V\7^ " (l + - l)- 1 ) W-\ri) < 0. (14) 

It remains to consider the case when = 0. It is easy to see that in this case the substitution 
^ = j] into the inequality (^1^1^) < £' < / l || x I'|| 2 leads to the condition even more stringent than 
Eq. (d. 

It makes sense to introduce the effective potential K//(y, R) = J <ix|</)o| 2 IV-. The function r\ 
has the factorized form r/ = </>o(x)/(y, R). From Eq. (|14j) we conclude that the system of four unit 
charges is unstable if for any fixed y 

pi 



1 + (V3/Wr) -I)" 1 K//>0. (15) 



We shall make one simplification, which helps carrying out all calculations analytically. We have 
W = Wi + W 2 , where W x = V u + V 24 and W 2 = V 13 + V 23 and obviously W_ < (PFi)_ + (W 2 )-. 
Breaking the kinetic energy term in Eq. (|15|) we deduce that the system is unstable if both of the 
following inequalities are satisfied p R — 4fi R ^1 + (a/3/ (8/7p) — l) -1 ^ > for i = 1,2, where 
= / dx|0o| 2 (VFi)_. Using the explicit calculation from H we get v}f f < (3/16)|R+ (1 - 6)y|" 2 
and v}f f < (3/16)|R - fey|~ 2 . It is known 10] that p| - X/R 2 > for A < 1/4. Thus both 



inequalities are satisfied if 3//r (l + (\/3/(8 / u/j) — l) _1 j < 1. Solving this inequality and rescaling 
the masses tells us that the system is unstable if fi R < (13 — 2\/22)M, r /54, which proves Eq. (^Q). 
One can improve the constant in Eq. (^Q) if one extracts everything from Eq. ()15|) . We preferred 
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to make the simplification by splitting W into two terms because this makes the whole derivation 
analytical. Let us also remark that Instability in Eq. means that there is no bound state neither 
below nor at the threshold. Indeed, if we would have H^f = E^h^f then, because one can choose 
Vl/ > in the ground state, we immediately get (ty\H\fy) < E t h\\^\\ 2 which was the starting point 
of our analysis. 
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